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Abstract 

Our chances to halt epidemic outbreaks rely on how accurately we represent the popu- 
lation structure underlying the disease spread. When analyzing global epidemics this force 
us to consider metapopulation models taking into account intra- and inter-community in- 
teractions. Recently Watts et al introduced a metapopulation model which accounts for 
several features observed in real outbreaks [Watts et al, PNAS 102, 11157 (2005)]. In this 
work I provide an analytical solution to this model, enhancing our understanding of the 
model and the epidemic outbreaks it represents. First, I demonstrate that depending on 
the intra-community expected outbreak size and the fraction of social bridges the epidemic 
outbreaks die out or there is a finite probability to observe a global epidemics. Second, I 
show that the global scenario is characterized by resurgent epidemics, their number increas- 
ing with increasing the intra-community average distance between individuals. Finally, I 
present empirical data for the AIDS epidemics supporting the model predictions. 

Human populations are structured in communities representing geographical locations and 
other factors leading to partial segregation. This population structure has a strong impact on the 
spreading patterns of infectious diseases among humans, forcing us to consider metapopulation 
models making an explicit distinction between the intra- inter-community interactions [Rvachev 
& Longini, 1985, Sattenspiel & Dietz, 1995]. The increase in model realism is paid, however, 
by an increase in model complexity. Detailed metapopulation models are difficult to build 
and as a consequence they are available for a few locations in the world [Rvachev & Longini, 
1985, Flahault et al., 1988, Eubank et al., 2004, Germann et al., 2006] or they cover a single 
route of global transmission [Hufnagel et al., 2004,Colizza et al., 2006]. 

Recently Watts et al [Watts et al., 2005] introduced a simple metapopulation model making 
an explicit distinction between the intra- and inter-community interactions. In spite of the model 
simplicity it accounts for several features observed in real epidemic outbreaks. In particular, the 
numerical results indicate the existence of a transition from local to global epidemics when the 
expected number of infected individuals changing community reaches one [Watts et al, 2005]. 
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I go a step forward and provide an analytical solution to the Watts et al metapopulation model. 
I demonstrate that there is indeed a phase transition when the expected number of infected 
individuals changing community reaches one. This analytical solution allow us to obtain a 
much deeper insight into the main features of global epidemic outbreaks. 

Model 

Figure [T] illustrates the general features of an epidemic outbreak on a population structured in 
different communities. Starting from an index case a disease spreads widely inside a com- 
munity thanks to the frequent intra-community interactions. In addition the disease is trans- 
mitted to other communities via individuals belonging to different communities. While the 
inter-community interactions may be rare they are determinant to understand the overall out- 
break progression. Based on this picture I divide the population in two types or classes. The 
locals belonging to a single community and the social bridges belonging to different commu- 
nities. In a first approximation I assume that (i) all communities are statistically equivalent, (ii) 
the mixing between the local and bridges is homogeneous, and (iii) social bridges belong to 
two populations. While these assumptions are off course approximations they allow us to gain 
insight into the problem. They could be relaxed in future works to include other factors such as 
degree correlations among interacting individuals [Vazquez, 2006c] and more realistic mixing 
patterns [Vazquez, 2006 d]. 

An epidemic outbreak taking place inside a community is then modeled by a a multi-type 
branching process [Mode, 1971] starting from an index case (see Fig. [l]). The key intra- 
community magnitudes are the reproductive number and the generation times [Anderson & 
May, 1991, Vazquez, 2006b]. The reproductive number is the average number of secondary 
cases generated by a primary case. The disease transmission introduces some biases towards 
individuals that interact more often. Therefore, I make an explicit distinction between the in- 
dex case and other primary cases and denote their expected reproductive numbers by R and 
R, respectively. The generation time r is the time elapse from the infection of a primary case 
and the infection of a secondary case. It is a random variable characterized by the generation 
time distribution function G(r). These magnitudes can be calculated for different models such 
as the susceptible infected recovered (SIR) model and they can be estimated from empirical 
data as well. Finally, a community outbreak is represented by a causal true rooted at the index 
case [Vazquez, 2006a, Vazquez, 2006b]. In this tree the generation of an infected case is given 
by the distance to the index case. Furthermore, the tree can have at most D generations, where 
D is the average distance between individuals inside a community. 

Spreading regimes 

Let us focus on a primary case at generation d and its secondary cases at the following gener- 
ation (see Fig. |2j). Let Nd(t) denote the expected number of descendants of the primary case 
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Figure 1 : Epidemic outbreak on a structured population: Schematic representation of a pop- 
ulation structured in communities (big circles) and and the spread of an infectious disease in- 
side and between communities. The individuals inside a community are divided between locals 
(open circles) and bridges (filled circles). The locals transmit the disease to other individuals 
inside the community (solid arcs) while the bridges transmit the disease to individuals in other 
communities (dashed arcs). For simplicity individuals that are not affected by the outbreak are 
not shown. 
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Figure 2: Local disease transmission: Diagram representing the disease transmission from a 
primary case at generation d to secondary cases in the following generation. The secondary 
cases are locals with probability 1 — (3, potentially leading to subsequent infections inside their 
community, and bridges with probability (3, transmitting the disease to other communities. Note 
that the expected number of descendants generated by a secondary case is evaluated at a delayed 
time t — T, where r is the generation time. 

at generation d. In particular N (t) gives the expected number of descendants from the index 
case, i.e. the expected outbreak size. In turn, N d+1 (t) is the expected number of descendants 
generated by a local secondary case at generation d + 1. Otherwise, if the secondary case is 
a bridge, it starts a new outbreak in a different community with expected outbreak size N (t). 
Putting together the contribution of locals and bridges we obtain the recursive equation 



(1-/9) 



1 + R f dG(r)N d+1 (t-r) 
Jo 



d = 



Nd ^ i (1-/3) [l +RJ* dG(r)N d+1 (t - r)l + (3N (t) , < d < D 



{ l-p + pN (t), 
Iterating this equation from d = D to d = we obtain 



d — D. 



,(t) = 1 + (1 - (3)F(t) + (3 [ dF(r)N (t - r) , 

Jo 



where 



(1) 



(2) 



D d _ l 

F(t) = Rj2[a-P)R] " G*(t) 

d=l 



(3) 



andG^(t) denotes the d-order convolution of G(t), i.e. G*°(t) = 1 and G* d+1 (t) = f* dG(r)G* d (t- 
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r). F(t) represents the expected outbreak size inside a community at time t and 



N c = lim F(t) , (4) 

t— »oo 

is the final expected outbreak size inside a community. When /3 = it coincides with the 
expected outbreak size inside a community [Vazquez, 2006&]. When /3 > © provides a self- 
consistent equation to determine the overall expected outbreak size after taking into account the 
inter-community transmissions. 

To calculate N (t) I use the Laplace transform method. Consider the incidence 

I(t) = N (t) (5) 

and its Laplace transform 

POO 

I(u) = / dte- wt I(t) . (6) 
Jo 

Substituting the recursive equation © in © I obtain 

ft \ /M m 
Iiuj) = , (7) 

where 

POO 

f{u) = / dte~^F{t) . (8) 



The validity of © is restricted to uo values satisfying 1 — (3f(u) > 0, resulting in different 
scenarios depending on the value of the parameter 

#c = PN C . (9) 

Local outbreaks: When Rc < 1 then I(u) is defined for all uo > and I(t) is obtained inverting 
the Laplace transform in ©. Furthermore, since 1(0) is defined from © it follows that I(t) 
decreases to zero when t — > oo, i.e. the epidemic outbreak dies out. 

Global outbreaks: When R c > 1 the incidence grows exponentially I(t) ~ e Wc< , where uj c is 
the positive root of the equation 

/?/>) = 1 . (10) 

These two scenarios are equivalent to those obtained for a single community [Anderson & 
May, 1991]. Rq represents the effective community's reproductive number and the threshold 
condition 

Rc = l (11) 
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delimits the local and global scenarios. 

To go beyond the final outbreak I analyze the progression of the inter-communities out- 
break. I assume that the disease is transmitted at a constant rate A from a primary case to a 
secondary case independently of their type. In this case the intra-community incidence is given 
by [Vazquez, 2006&] 



X(Xt) 



F(t) ~ N C — , (12) 



for t ^> to, where 

D — 11 



R X 

Calculating the inverse Laplace transform of © I finally obtain 



(13) 



oc 



\( \ f\D(m+l)-l p -\t 

m = N„j:w-r ( r [' g(m+1)1 ■ 

where T(x) is the gamma function. Figure|3]shows the progression of the incidence as obtained 
from (fT4l) . As predicted above, the outbreak dies out when R c < 1 while when R c > 1 it 
grows exponentially. More important, the incidence exhibits oscillations at the early stages, 
their number increasing with increasing D. For example, we distinguish about two oscillations 
for D = 10 while for D = 30 several oscillations are observed. These oscillations represent 
resurgent epidemics, which are often observed in real outbreaks [Riley & et al, 2003, Anderson 
& et al, 2004] and simulations [Sattenspiel & Dietz, 1995, Watts et al, 2005]. 



Case study: AIDS epidemics 

To understand the relevance of these results in a real world scenario I analyze data reported 
for the AIDS epidemics. First, I estimate the parameter Rc determining the spreading regime, 
local or global. Figure |4^ shows the value of Rc across the USA by state. For most states 
Rc > 1, reaching significantly large values for several states. For example, Rc exceeds 1,000 
for California and New York. These numbers indicate that the USA AIDS epidemics is in the 
global spread scenario (Rc > 1), in agreement with the general believe. 

Second, I analyze the temporal evolution of the AIDS incidence. Fig. and c show the 
AIDS incidence in USA and UK by year, indicating a similar temporal pattern. The epidemics 
started with an increasing tendency of the incidence which, after reaching a maximum, switched 
to a decreasing trend. After some years, however, the epidemics resurges with a new incidence 
increase. This picture coincides with the model predictions in Fig. |3] Therefore, a possible 
explanation of the observed multiple peaks is the existence of a community structure, which 
can be attributed to geographical location and other factors. 
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Figure 3: Epidemic outbreak progression: The incidence I(t) as a function of time in units 
of the local disease transmission rate A, for R c = 0.1 (dash-dotted), 0.5 (dotted), 0.9 (dashed), 
1.0 (solid) and 1.1 (dash-dash-dotted). The panels from top to bottom corresponds to different 
average distances D between individuals inside a community. 
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Figure 4: USA AIDS epidemics: Estimated R c = (3N C for the different USA states. f3 was 
computed as the ratio between the number of state out-immigrants and the total state population 
according to the 1995-2000 USA census (http://www.census.gov). N c was computed as the 
number of habitants living with AIDS according to the 2005 statistics published by the US 
Department of Health (http://www.hhs.gov). b) and c) AIDS incidence in the USA b) and UK 
c) by year, as reported by the US Department of Health and the UK Health Protection Agency 
(http://www.hpa.org.uk), respectively. 
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Discussion and conclusions 



Rc in (ITU represents the expected number of infected individuals leaving their community. 
The numerical simulations reported in [Watts et ah, 2005] indicated the existence of a transi- 
tion at R c = 1, from local outbreaks when R c < 1 to global epidemics when R c > 1. I 
have demonstrated that there is indeed a phase transition at R c = 1. Furthermore, the an- 
alytical solution provides an expression of Rc as a function of the bridge's fraction and the 
intra-community expected outbreak size ©. Rq represents a measure of the reproductive num- 
ber at the inter-community level. Its value can be estimated from the expected outbreak size 
inside a community and the bridge's fraction. Based on the resulting estimate we can determine 
if an epidemics is in the local or global epidemics scenario and react accordingly. 

The inter-community disease transmission is characterized by oscillations at the early stages 
which represents resurgent epidemics, the number of these resurgencies being determined by 
the characteristic distance between individuals within a community. In essence, when D is 
small the time scale characterizing the outbreak progression within a community is very small 
[Barthelemy et ah, 2004, Barthelemy et ah, 2005, Vazquez, 2006Z?]. Therefore, the time it takes 
to observe the infection of a social bridge is very small as well, resulting in the mixing between 
the intra- and inter-community transmissions. In contrast, when D is large it takes a longer time 
to observe the infection of a social bridge and by that time the intra-community outbreak has 
significantly developed. Therefore, in this last case the outbreak within communities is partially 
segregated in time. 

When multi-agent models are not available these results allow us to evaluate the potential 
progression of an epidemic outbreak and consequently determine the magnitude of our response 
to halt it. They are also valuable when a detailed metapolpulation model is available, funneling 
the search for key quantities among the several model parameters. More important, this work 
open avenues for future analytical works that side by side with multi-agent models will increase 
our chances to control global epidemics. 
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